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Abstract 



We study the Fermi gas quantum mechanics associated to the ABJM matrix model. 
We develop the method to compute the grand partition function of the ABJM theory, 
and compute exactly the partition function Z{N) up to = 9 when the Chern-Simons 
level k = 1. We find that the eigenvalue problem of this quantum mechanical system is 
reduced to the diagonalization of a certain Hankel matrix. In reducing the number of 
integrations by commuting coordinates and momenta, we find an exact relation concern- 
ing the grand partition function, which is interesting on its own right and very helpful 
for determining the partition function. We also study the TBA-type integral equations 
that allow us to compute the grand partition function numerically. Surprisingly, all of 
our exact results of the partition functions are written in terms of polynomials of tt~^ 
with rational coefficients. 



July 2012 

*yhatsuda(§th . phys . titech .ac.jp 
^mor iy amaSmath . nagoya-u .ac.jp 
"'"kazumiSazusa. shinshu-u. ac . jp 



1 Introduction 



M-theory is expected to unify all of our understandings on the nonperturbative effects of 
string theory. Despite its importance, many aspects on this interesting theory await to be 
understood. One of the most mysterious facts on M-theory is the degrees of freedom of the 
branes. From the AdS / CFT correspondence, it is known [Tj that coincident M2-branes (or 
M5-branes) have the degrees of freedom of order N^/"^ (or 7V^, respectively). In the case of 
D-branes in string theory, the worldvolume theory can be described by the super Yang-Mills 
theory, which has the degrees of freedom A^^. Compared with such a case, the field theoretical 
realization on the branes in M-theory is still unclear. 

Recently, there is a breakthrough [2] (along with its extension with fractional branes [3]) in 
understanding branes in M-theory, where it was proposed that the worldvolume theory of M2- 
branes is described by the supersymmetric extension of the Chern-Simons theory. Concretely 
speaking, the proposal states that A^ coincident M2-branes on the geometry M'^/Z^ is described 
by the M = Q Chern-Simons theory with gauge groups U{N)k x U{N)^k and bifundamental 
matters Ai^2, Bi^2 forming the superpotential W = ^ Tt{AiBiA2B2 — A1B2A2B1) where the 
coupling constant is given by = ^. The geometry R^/Zk in a special case with /c = 1 is 
nothing but the flat space, hence it is natural to expect that the supersymmetry is enhanced 
to A/" = 8 in that case. It was then proposed that assuming the existence of the monopole 
operator we can construct the A/" = 8 supersymmetry algebra when /c = l,2[ll|5l|6l[7]. 

Subsequently, it was shown in [8] that the localization technique [9], which reduces the 
partition function defined by infinite-dimensional path integral to a finite-dimensional matrix 
model, is applicable to supersymmetric Chern-Simons theories. This matrix model has the 
measure of the Chern-Simons type and many techniques from the Chern-Simons matrix model 
such as [TU] was utilized [TT] to obtain not only the N^/"^ behavior of the free energy in the 
planar limit but also the shift of 't Hooft couphng x T ~ ^ expected from the Euler 
coupling in the gravity dual [T^ . 

Besides, it was also noticed that this matrix model is closely related to the topological string 
theory on a local Calabi-Yau over the Hirzebruch surface Fq = x P^. Using this relation, 
we can study the higher genus contribution by borrowing a celebrated recursive relation, the 
holomorphic anomaly equation [13], from the topological string theory. Following interesting 
observations by [13], it was found in [15] that we obtain the Airy function after summing up 
the contributions from all genus. 

This result on the Airy function suggests some interesting relations to previous under- 
standing on M-theory. In [16] it was noticed that the Airy function also appears in other 
context of M-theory as the wave-function of the universe. Further, the integral expression of 



1 



the Airy function Ai (z) = J^dtexp{—zt + y), which is a contour integral of an exponenti- 
ated cubic term, looks similar to the partition function Z = J DAexp{—S[A]) of the cubic 
Chern-Simons action S[A]. Therefore, we are led to expect a deep relation among M-theory, 
the Chern-Simons theory and the Airy function. 

Another expectation from this result is its implication in understanding some of hidden 
structures in the membrane theory. String theory originates from the Veneziano amplitude 
in the dual resonance model, out of which the string worldsheet conformal symmetry was 
realized. Similarly, one might expect to find some hidden symmetry out of this Airy function 
and realize some structure for M-theory. An interesting hidden structure of this Airy function 
was beautifully realized by [17]. It was shown that the partition function of M2-branes (as well 
as its generalization with = 3 supersymmetry) can be rewritten exactly into the partition 
function of a Fermi gas system with the density matrix (related to Hamiltonian by p = e~^) 



^2 cosh I 2 cosh I ^2 cosh f ^'^ 

and, with the standard analysis of the statistical mechanics, we can reproduce the behavior 
of N^/'^ and Airy function easily. Very recently, it was further proposed in [18] that we can 
generalize this formalism to the M = 2 case by introducing interactions to the Fermi particles 
and in [19] that we can extend to the cases with Wilson loops. 

This simple rederivation of the Airy function implies the importance of this Fermi gas for- 
mulation. Since we have already a good control on the perturbative terms using the statistical 
mechanical method, as a next step, we would like to proceed to the study of this quantum 
mechanics extensively beyond the perturbation theory and to obtain some non-perturbative 
information out of it. 

Aside from its possible important message on M-theory, the mathematical structure in 
this maximally supersymmetric theory seems very interesting. The supergravity analysis [20] 
suggests the BPS condition is expressed by a continuous Toda equation and the Hamiltonian of 
this Fermi gas system takes the form similar to the relativistic Toda theory [211 122] . Therefore, 
it is natural to expect some kind of the exactly solvable structure in this theory. Some of the 
exact results was found in [23] . 

We would like to make a first step towards the goal to compute the exact partition function 
of the ABJM theory for arbitrary and arbitrary Chern-Simons level fc, and to solve the 
eigenvalue equation of the corresponding Fermi gas system. Heading for this goal, we have 
found some interesting exact results on this theory as we will report in this paper. 

As we will see in the next section, the (grand) partition function of this Fermi gas system 
is computed by the trace of the power of the density matrix. The main difficulty resides in 



2 



the computation of convoluting the matrix element successively, 

p''{Qo,Qn)= /^^■■■%^p(go,gi)p(gi,g2)---p(gn-i,gn). (1-2) 

J Zn Zn Zn 

Following the previous work on the thermodynamic Bethe ansatz (TBA) equations [211 ESj , 
here we are able to express the matrix elements in terms of more basic vector quantities and 
reducing the computation of convoluting matrices into that of convoluting vectors. 

In addition, we have also attempted to reduce the integration by commuting the coordinate 
q and the momentum p. We have found that, although the naive commutation relation 
is anomalous, interestingly the additional term to the naive commutation relation can be 
expressed in terms of a projection operator \p = 0){q = 0\. Using this relation, we can express 
Trp" for odd n by some combinations of {q = 0\{2coshq/2)~^^'^p"^{2coshq/2Y^'^\p = 0). For 
even n, we have found an interesting novel relation which separates the density matrix into 
two parities. 

We will mostly restrict ourselves to the maximally supersymmetric case k = 1. Most of 
our following analysis goes straightforwardly to other values of k though a key relation will 
only look simple for k = 1. The case k = 1 is particularly interesting from the viewpoint of 
holography since the ABJM theory at /c = 1 is conjectured to be dual to the M-theory on 
AdS^ X S*^ which is in the truly M-theoretic regime. The study of the ABJM matrix model in 
the large A^, finite k regime was initiated by [SB] and further studied in [T7| from the Fermi 
gas approach. However, the approach of [T^ is restricted to the perturbation around k = 0. 
We would like to study the k = 1 case exactly. 

As a result, we have computed the exact partition function Z{N) up to A^ = 9. We shall 
summarize our results in the following, 

Z{1) = I Z{2) ^ 



4' ^ ^ 167r' 

. s 7r-3 10 -vr^ 

-97r2 + 207r + 26 , SQtt'^ - Uln^ + 7^ 



-7571^ + IQStt^ + 1747r - 126 _ IOSStt^ - 20167r3 - 41487r2 + 876 

SSITtt^ - 134807r3 - 153487r2 + 88807r + 4140 
^(9) • ^^■^> 

Surprisingly, all of these results are written in terms of polynomials of vr"^ with rational 
coefficients. 

This paper is organized as follows. In Section 2, we present our general procedure to 
compute the grand partition function. Applying the lemma of Tracy- Widom [24], we first 
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show that the power of density matrix q') is written in terms of a set of functions of one 
variable. We also show that the eigenvalue problem of p is reduced to the diagonalization of a 
certain Hankel matrix. In Section 3, we find that the commutation relation of (2 cosh |)~^ and 
(2 cosh |)~^ appearing in p contains an anomalous term, which can be expressed by a insertion 
of projection operator. Then we show that such anomalous commutation relation leads to a 
nontrivial relation of the grand partition function. In Section 4, using our general method, we 
compute exactly the first few terms of the Taylor expansion of the grand partition function 
at = 1. These are nothing but the partition functions Z{N). Along the way, we find a 
non-trivial equation relating the even part and odd part of density matrix. In Section 5, some 
numerical results are given. We consider the TBA-type integral equation proposed in [25] to 
compute Tr p". These integral equations are powerful in the numerical analysis. We also study 
the numerical solutions of the eigenvalue problem of the Fermi gas system. In Section 6, we 
compare our results with the known perturbative results of the ABJM matrix model and find 
a good agreement with the Airy function even for small A^. We study the non-perturbative 
corrections to the free energy, and find numerical evidence that the leading correction scales 
as e"^'^^'^. We also consider the expectation value of A^ in the grand canonical ensemble. We 
find that the numerical computation using TBA-like equations smoothly interpolates our exact 
results and the perturbative results. Finally, we conclude in Section 7 with discussions and 
future directions. Appendices A, B and C contain some details of the computations skipped 
in the text. 



2 Grand partition function of ABJM theory 

The aim of this section is to present a general prescription which is helpful for computing the 
grand partition function of the ABJM theory. The concrete computation based on the method 
given here will be presented in Section |H 

Let us start by reviewing the Fermi gas formalism developed in [T7]. By using the local- 
ization, the partition function of the ABJM theory on is written as the matrix integral 



^ J_ f d^p d^y n.<,[2sinh(^)]^[2sinh(^)]^ 
^> J {27r)^ {27t)^ n,,[2cosh(^)]2 ''''^ 



4 = 1 



(2.1) 

Applying the Cauchy identity to the measure factor and performing the Gaussian integral over 
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fii, z/j, this partition function is recast into the "mirror" expression f27\ [IZ] 



ZiN) 



2^N\ J (27rA;)^ J-i 2cosh 



n 9 rnsh 21 11 



tanh 



i<j 



2 Q i - Qj 
2k 



The key idea of [T7] is to rewrite it as the partition function of an ideal Fermi gas 



■ o-GSjv ^ ' i=l 



where 



(2coshf)V2 2cosh(2i^) (2coshf)V2^ 



(2.2) 



(2.3) 



(2.4) 



is understood as the density matrix of this Fermi gas system. Since the sum over permutations 
in fl2.3p can be also written as a sum over conjugacy classes of the permutation group, it is 
convenient to consider the grand partition function 



1 + Z{N): 



.N 



(2.5) 



N=l 



where 2; = e'^ is the fugacity and /i is the chemical potential. In our analysis, it is important 
to notice that the grand partition function has the following form 



where Tr is defined by 
Trp" 



[z) = exp 



dqi dqr, 



E 

n=l 



n 



Trp" 



(2.6) 



27ik 27ik 



p{Qi, Q2)p{q2, qs) ■ ■ ■ piQn, qi)- 



(2.7) 



It is easy to see that ^(2) can be expressed as the Fredholm determinant of the kernel p0 



E{z) = det(l + zp). 



{21 



Once the grand partition function is known, one can obtain the partition function by perform- 
ing the contour integral 



Z{N) 



dz 2(2) 
2-Ki 2^+1 ■ 



(2.9) 



similar analysis of the grand partition function of other matrix models has been appeared in [28 1 129 1 130 ] . 
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For later analysis, it is convenient to divide the density matrix into two parts. Since the 
density matrix p{qi, ^2) is parity-preserving: p{—qi, — ^2) = pili, (I2), we decompose it into the 
parity even and odd parts 

p{(li,'l2) = P+{qi,q2) + P~{qi,q2), P±{qi,q2) = ^ • (2.10) 

By definition, p±{qi, — ^2) = ^P±{qi, 92), and thus all the matrix elements of the product p+p_ 
vanish: 

/"°° dq 

(p+P-) {(11,(12) = J — p+(gi,g)p_(g,g2) = 0. (2.11) 
This means that Tr p" splits into two parts, 

Trp" = Trp!^ + Trp!^. (2.12) 
and the grand partition function is factorized into 

S(z)=S+(^)S_(^), (2.13) 

where 

S+(z) = det(l + 2p+), S_(z) = det(l + zp_). (2.14) 

As we will see later, such a decomposition of the density matrix plays an important role in 
the analysis of the grand partition function. 

2.1 Method to compute the grand partition function 

Now let us consider the computation of the grand partition function. As can be seen from 
(12.61) . we can compute the grand partition function if we know Tr p". To compute the partition 
function Z{N), in particular, we need Trp" (1 < < A^). However it is technically difficult 
to directly compute Tr p" (or the matrix element p"(q'i, q'2)) because it is given by the multi- 
variable integral as in (12. 7p (or (II. 2p ). To overcome this difficulty, we here use the technique 
proposed in [21], which enables us to compute p"'{qi,q2) from a series of functions of one 
variable. This technique can be applicable to a wide class of the kernels with the form 

E{q,)E{q2) 

K{qi,q2) = ^ , y 2.15 

M(gi) + M{q2) 

Here we apply it to the kernels p±. Below we concentrate our attention to k = 1 for simplicity. 
It is straightforward to generalize the method here to general k >2. 
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The important point is that two matrices p±{qi,q2) have the following form; 

E±{qi)E^iq2) 



P± (91,^2) 



cosh qi + cosh q2 



where 



^ ^ ^ coshf ^ ^ ^ sinhf 



2 cosh I 



(2.16) 



(2.17) 



We can apply Lemma 1 in [21] to (I2.16p . It is shown that (the generating functions of) the 
density matrices are given by 

Q±iqi)Q±{q2) - P±iqi)P±{q2) 



zp± 
1 - z^pI 

z'pI 

1 - z^pI 



(91,92) 
(91,92) 



cosh qi + cosh q2 

Q±{qi)P±{q2) - P±{qi)Q±{q2) 
cosh qi — cosh q2 



where 



[dq' 


zp± 


1 2n 


_i-z^pI_ 



,q')^zE^{q'), Q^{q) = 



[dq' 


1 


1 2n 


i-^VL 




00 



(2.18) 
(2.19) 

,q')^E^iq'). 

(2.20) 



If we expand P±{q) and Q±{q) as 

00 

P±(g) = ^zE^{q)Y.z''^^'ct>'t\q), Q^{q) = v^i?±(9) ^''4'(9), 

where 



(2.21) 



k=0 



k=0 



4(9) 



dQ' 



pUq,q)Ei{q'), 



E±iq)J 27r' 

then (I2.18P and (12.19^ indicate that the density matrices are given by 

E±{qi)E^{q2) 



(2.22) 



p|"+'(9i,92) 
P± (91,92) 



cosh qi + cosh q2 

E±{qi)E^{q2) 
cosh qi — cosh q2 



E(-i)'>i(^i)'^"fe), 



fc=0 

271-1 



5^(-l)Vi(9i)0l"-'-'(92). 



(2.23) 
(2.24) 



k=0 



^Note that 



^(91,92) = 



itself is also of this class. 



e 



e 2 ■ 



^2 cosh § e9i + e92 ^2 cosh f ^2 cosh ^ e-9i + e-92 ^2 cosh f 
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These are the key relations in our analysis later. These results show that the matrix elements 
P±('?i) ^2) are determined algebraically by the functions (/c = 0, . . . , n — 1) of one variable 

respectively. Therefore our problem reduces to finding (p^^lq), which is much simpler than the 
original one of finding p^{qi,q2) themselves. 

For the determination of as discussed in Appendix|Al we can show that 4>'±{q) satisfy 

the following recurrence relations 



<Pliq) = c, I ^c,.,<p'-\q'), <pt{q) = c,j^^ %c,.,t\,<\>\^\^). (2.25) 

Throughout this paper, we often use the notations 

1 1 g 

Co = — r, Sn = --^7, ta = tanh -. (2.26) 

^ 2 cosh I' ^ 2sinh|' '2 ^ ' 

Prom (12.251) with the initial condition (p^{q) = 1, one can recursively fix 0j.(q'). In Section HJ 
we will indeed determine the first few 4>±{q) by using these recurrence relations. 

Let us summarize our general prescription to compute the grand partition function. We 
first determine (p^iq) by using (I2.25p . Then from (12.231) and (I2.24p . we find the matrix elements 
p±{qi,q2)- Computing Trp![. from such expressions, we obtain the grand partition function 
through ^M) and fl212D . 



2.2 Eigenvalue problem and Hankel matrix representation 

In the previous subsection, we considered Tr p" to know the grand partition function. There 
is another approach to analyze it. The grand partition function can be expressed in terms of 
the eigenvalues of the density matrix p, 

00 

E{z) = det{l + zp) = Y[{1 + zXr,), (2.27) 

n=0 

where A„ is the n-th eigenvalue of p with decreasing order Aq > Ai > ■ ■ ■ . Thus if we can 
solve the eigenvalue problem of p, we obtain '^{z). The eigenvalue equation is given by the 
homogeneous Fredholm integral equation of the second kind 

J ^P{q,q')^{q') = ^^{q). (2.28) 

It is, of course, hard to solve this equation, but we can rewrite it as the eigenvalue equation 
for an infinite dimensional Hankel matrix M whose matrix elements are given by 

1 /•°° tanh"+'"(^) 1 r+'" 
Mnm = — r / dq = — / dt , 2.29 

87Tkj_^ \osh(|)cosh2(4) 4ttJ_, mi/vT^y ^ ' 
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where Tk{x) is the k-th Chebyshev polynomial of the first kind. This representation is useful 
in our later analysis. 

To rewrite the eigenvalue equation, let us first expand the kernel fl2.4p into 



n=0 



where 



Pn(g) 



tanh'^(^) 



(2.30) 



(2.31) 



2^A; cosh(f)cosh(4)' 

Since Pn{—Q) = {^^)"Pn{Q), the decomposition ( l2.1Up indicates that p± are respectively ex- 
pressed as 



n=0 



P+(gi,g2) =^p2n{qi)p2n{q2), P-(gi,g2) = p2n+l{qi) p2n+l{q2 

n=0 

Plugging (rai]) into we obtain 

oo 

^fnPn(g) = A^(g), 
n=0 

where 

^Pn(g)^(g). 



(2.32) 



Substituting f l233D back to f l234D . we find 

oo 



(2.33) 



(2.34) 



(2.35) 



m=0 



where the matrix elements of M are given by (12.291) . The equation (I2.35P formally shows that 
all the eigenvalues of p are also the eigenvalues of the infinite dimensional matrix M. 

The matrix elements Mnm depend only on + m, and such matrices are called Hankel 
matrices. For the present case, it is obvious that Mnm = if n + m is odd. Namely, this 
Hankel matrix M has the following form 

...\ 



M 








717,1 





7712 








mi 





m2 





ma 


nil 





7712 





7773 








7712 





77l3 





m4 


m2 





7713 





m4 


















(2.36) 
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We shall call such a matrix a parity-preserving Hankel matrix. As can be seen easily, the 
parity-preserving Hankel matrix can be decomposed into two blocks of Hankel matrices with 
opposite parity 





nil 


1712 


mi 


m2 




m2 







...\ 



v 



M_ 



7 



nil 


ni2 


m-s 


ni2 




m^ 




m4 





A 



(2.37) 



These two Hankel matrices are related by 

M_ = M.T^ 



(2.38) 



with the shift matrix {T±)mn = Sm±i,n- The eigen-spaces of M also decompose into the direct 
product of the eigen-spaces of M±. 

For k = 1, Mnm has the following very simple form {n + m: even) 



where C„ is the Catalan number 



2 

2»i+m,+3 ' 



{2n)\ 



[n + lV.nl 



(2.39) 



(2.40) 



For k = 2, we find (n + m: even) 



8tt 



n + m + 1 



n + m + 3 



n + m + 1 



(2.41) 



where il){x) = r'(x)/r(x) is the digamma function. 

After seen the simple Hankel matrix representation, before closing this subsection, now 
let us express various quantities appearing in the previous subsection in terms of this Hankel 
matrix M. Since the matrix M and the density matrix p share the same set of eigenvalues 
assuming a certain regularity condition, it is obvious that Tr is equal to Tr p" for any n, 
thus the grand partition function is given by 



E{z) = det(l + zM). 



(2.42) 



Repeating the same argument above, one finds that the eigenvalue problems for p± are mapped 
to those for the Hankel matrices M±, respectively. Let A-t „ {n = 0,1,2,...) be the n-th 
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eigenvalue of the density matrix p±. Then, we find 

oo 

JJ(l + 2A+,„) = det(l + 2:M+), 

n=0 

oo 

+ 2A„,„) = det(l + 2M_). (2.43) 



oo 

s+W = 

n=0 

oo 

n=0 



We observe from the numerical evaluations (see Section 15.31) that the eigenvalues of p± satisfy 
the inequality 

> A_,o > > A_,i > A+,2 > A_,2 > ■ ■ • , (2.44) 
which means the following relations 

A+,n = A2r!,5 A_^„ = A2n+1- (2.45) 

Also, (l)\{q) and its Fourier transformation (j)\{p) 



(p) = J ^e-^Pyliq). (2.46) 
can be expressed in terms of M as 

oo oo 

lC^"<P'M= EC(^+)-' 2c-^<^t(g) = EC(^+"')"0. (2.47) 

m,n=0 n=0 

where tq = tanh(g/2). The derivation of these representations will be given in Appendix [Bl 
2.3 Representation in Chebyshev basis 

In the previous subsection, we have found an interesting Hankel matrix representation for 
the density matrix p{qi,q2)- It is natural to wonder whether there is an orthonormal basis 
in which the density matrix p{qi,q2) is represented by this Hankel matrix. We will show in 
this subsection that the Hankel matrix M is closely related to the density matrix in a basis 
constructed by the Chebyshev polynomials. 

Motivated by the expansion (12.301) . let us introduce the following function 

/ I \ V2^n(tanh(|)) 

\^\^) = , ' (2-48) 



cosh3/^( 



where Un{x) is the Chebyshev polynomial of the second kind and its appearance originates 
from the orthonomality of the basis {\n)} 



oo 

n=0' 



dq v/2y^(tanh(|)) v/2[/„.(tanh(|)) _ 
^ ~ ' cosh=^/^(|) cosh^/^d) " ^ ^ ^ 
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Then, the density matrix in this basis, Mp_^^^ = {ni\p\n2), is given by 

^nma = / ^^{ni\qi){qi\p\q2){q2\n2) = ^ / ^ {qi\ni) pn{qi) j ^(g2|ri2)Pn(g2) 

n=0 



(2.50) 



where we have used the expansion ( ]2.30p . Let B be the matrix for the basis change from 
Chebyshev polynomials to monomials 

oo 

r = Y,BnmUUt)- (2.51) 

m=0 

The explicit form of the non- vanishing matrix element Bnm is given by 

B„„= l!!i±lf, " + 1 J. (2^52) 



2" n + 1 V('^ - 'm)/2 
when ?T, + m is even and n> m. Since one can check 

the matrix element M^^^^ is given by 

1 °° 1 
K.n. = g E ^-1^-2 = -{B^BUn,. (2.54) 

n=0 

On the other hand, the Hankel matrix M for /c = 1 is written as 

1 f tanh"i+"^(|) 1 o o f . f/mi(tanh(|))f/„,(tanh(|)) 



1 /• tanh"i^"2(f) _ 1 ^ [ 
- ^ y cosh=^(£) ~ 8^ ^ B^,m,B^,^, j c 



. cosh=^(f) 

■•2' mi,m,2 ^2' 

= ^(i?5^)nin.. (2.55) 

In summary, the Hankel matrix M and the density matrix in the Chebyshev basis are 
written as 

M = -BB^, MP = -B^B. (2.56) 
8 8 

These two matrices are isospectral, i.e., they have the same set of eigenvalues (assuming a 
certain regularity condition) 

Mv = Xv ^ M^v' = Xv', v' = B^v. (2.57) 

Note that is not a Hankel matrix while M is so. 
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3 Algebraic treatment 

Here we consider an algebraic treatment of our Fermi gas system. Note that the matrix 
element in the position representation (12 ■4p 

p(?i,92) = (gi|p|g2), (3.1) 

comes from a one- dimensional quantum mechanical system with the operator p expressed in 
terms of the coordinate q and the momentum p as 

p = y/c^Cp^. (3.2) 

Below we focus on the k = 1 case. Let us summarize our convention for the phase space used 
in this paper, which is especially clean for k = 1. The commutation relation of q and p is 
given by 

[q,p] = 27n. (3.3) 
The normalization of the bases \q) and \p) is as follows: 

{q\q') = 27r6{q - g'), {p\p') = 2vr5(p -p'), (qIp) = e^^'. (3.4) 
Therefore the projectors are normalized as 

Hereafter, we suppress the hats on operators for simplicity. Note that the matrix elements of 
Cp and Sp in the position representation are 

{qi\cp\q2) = Cq,-g^, (gi|sp|g2) = ^tanh (^^^-^^ • (3-6) 
Using the bra-ket notation, the matrix elements of p" are given by 

P"(gi,g2) = {qi\{V^,Cp^r\q2) = .[^{qi\{cpcX\q2). (3.7) 



Since Tr p" = Tr[(cpCg)"] = Tr[(cqCp)"] from this equation, the grand partition function is also 
written as 

'^{z) = det(l + zCpCq) = det(l + zCqCp). (3.8) 
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3.1 Anomalous terms in commutation relations 

As explained around fll.2p . the difficulty in computing Tr p" consists in the multi-dimensional 
integral. In the operator representation, it can also be stated as the difficulty in the integration 
over the multi-dimensional phase space. Therefore, we would like to reduce the dimension. 

If we know the commutation relation between Cp and Cq, we can reduce the dimension of 
the phase space involved. As is seen here, one should be careful to derive the commutation 
relation. From the operator identity 

2 cosh - ■ 2 cosh - = i ■ 2 sinh - ■ 2 sinh -, (3.9) 

2 2 2 2' ^ ^ 

one would naively expect 

2cosh| ■ 2coshf ~ 2sinh| ' 2sinhf ^^'^^^ 

However this is not the case. There is an anomalous term in this commutation relation. To 
see this, let us consider the matrix element of the difference between both sides of (13.101) 

1 1 1 1 1 q2 



If we introduce a projection operator \p = 0){q = 0| =: |0p)(0g| satisfying (|0p)(0q|)^ = |0p)(0 
we can express the relation as 



CpCg = -isgllsp, n = 1 - |0p)(0q|, (3.12) 

at the operator level. Similarly, we find 

CgCp = iSpUUg, nt = 1 - |0,)(0p|. (3.13) 

Using the commutation relations fl3.12p and fl3.13p . we can calculate Tr p" correctly. Let 
us see the computation of Tr p^ as a simple example: 

Trp^ = Tr^CpCgCpCg) = -iTi^SgllspCpCg) = -iTT{s2gIls2p), (3.14) 

where we have used SgCg = S2g. Since 11 = 1 — |0p)(0g|, (13.141) is further rewritten as 

Trp^ = J ^^S2gS2p + i J ^{q\s2g\0p){0g\s2p\q). (3.15) 

The first term vanishes because the integrand is an odd function. The second term is evaluated 
as 

/dq 1 f dq 2 

^{q\s2qK){0g\s2p\q} = -j (3.16) 
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where {q\s2q\0p) = S2q and {0g\s2p\q) = — (z/4) tanh(g/4). Therefore we obtain 

Tr/ = ^. (3.17) 

One can easily check that the same result is reproduced by a brute force integration using 
(12.71) . In the above computation, it is obvious that the anomalous term in the commutation 
relation is crucial to obtain the correct answer. In principle, we can compute Tr p" by using the 
relations (13.12^ and ( 13.13^ repeatedly instead of doing multi-integral. In practice, however, the 
computation becomes very complicated as n grows. One advantage to consider the algebraic 
treatment is to give us some non-trivial implications for the grand partition function as we 
see in the next subsection. 



3.2 Implications for the grand partition function 

In the previous subsection, we originally intended to reduce the dimensions of the phase 
space integration by commuting the coordinate q and the momentum p. Though this was 
not successful because of an anomaly term, the anomaly term is controlable and gives us 
an algebraic method to compute Tr p^. In this subsection, we would like to discuss a more 
important implication. 

Using (I3.13P and (13.121) . it is not difficult to show 

Tr(cgCpn)" = TrispIlUgU)'' = i"" TiiliUqIispY = i^"" TiiU^CpcX- (3-18) 

Since Tr(cqCpn)" is real, this should be equal to its conjugate TT{WcpCq)"'. The observation 
(I3.18P shows that, if the power n is odd, we have Tr(cgCpn)" = which implies that, by ex- 
panding n = 1 — I Op) {Oq\, Tr p" = TT{cqCp)^ can be expressed in terms of a certain combination 
of {Oq\{cqCp)"^\Op). The easiest way to incorporate the combinatorics correctly is as follows. 
Let us consider the determinant Si(z) = det(l + zCpUcq). This determinant is rewritten as 

E,{z) = det (l + zcpcq - ^|0p)(0,|) , (3.19) 

where we have used Cq\Oq) = (1/2) |0g). This determinant is formally written as 

Ei{z) = det(l + zCpCq) det ( 1 - ^- ^ |0p)(0,| ) . (3.20) 

V 41+ ZCpCq J 

The first factor is just the grand partition function S(2;), and (the logarithm of) the second 
factor is evaluated as 

(1 \ °° 1 / 1 \" 

^-iTT7jrj;l"'><«.lj=-T''E;r(iTT^W(".l) • P.^D 
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Due to the insertion of the projection operator, the last equation is drastically simplified 

1 



1 / 1 

n=l ^ ^ ^ 

Therefore the relation (13.201) becomes 

where we have introduced the new function G{z) by 



Op) 



(3.22) 



(3.23) 



G(z) = l + -(0, 



1 ZCpCq 



1 ZCqCp 



I Op). 



(3.24) 

Let us show that Hi (2;) is invariant under z — )■ —z. Using the commutation relations (I3.13P 
and (I3.12p . we have 

Ei{z) = det(l + zUcgCp) = det(l + zUiSpW^Sg) = det(l - zCpCgW^). (3.25) 
Since Hi (2) is a real function for 2 G M, this implies 

Hi(z) = det(l - z{cpCgU^y) = det(l - zUcgCp) = Ei{-z). (3.26) 
Plugging (13.231) into (13.261) . we find a relation 

^{-z) G{-z) 

Now note that the left hand side is nothing but (the square of) the grand partition function 
with odd power terms. More explicitly, let us divide ^(2) into two parts 



:i{z) — ^odd(2;)^even(^) 



n=l 



,2n-l 



2n - 1 



where 

Hodd(2) = exp 
By construction, one easily finds 

Then, (I3.27P shows 



Tr p 



2n-l 



'^even(^) — CXp 



n=l 



2n 



Tr p 



2n 



(3.28) 



(3.29) 



'='odd(2;) 



-z 



'='odd(2j 



m _ G{z) 
E{-z) G{-zy 



(3.30) 



(3.3i; 
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Using this relation, we can compute Tr p^""*"^ from G{z). As we will see in the next section, 
the function G{z) also has an important meaning for the even power terms, which connects 
the two determinants '^+{z) and defined in the previous subsection. 

After introducing the bra-kets such as (Op|, we can simplify various definitions and relations 
appearing previously. Using the bra-ket notations, the functions and (j)\{q) are expressed 

as 

= {q\M^\%). 4>1iq) = {q\M'\Og). (3.32) 
It is easily shown that (p^iq) satisfies the recurrence relation fl2.25p 

= / %{QMq'){q'\M''-'K) = CgJ^c,_g,<j>1-\q'). (3.33) 

Similarly, 0+(g) also satisfies the recurrence relation 



''{1) = j'%cg.g,cg4l'\q'). (3.34) 

Finally let us comment on some useful algebraic relations steming from (12.231) and (12.241) by 
setting one of the coordinate to zero. First, comparing (I3.32p with (13. 7p . one finds the relation 
between p'' and 0^ 

p''{q,0) = ^g4>1{q)- (3.35) 

Since p^(g,0) = for A; > 1, p^{q,0) = p+(g,0) holds. Secondly, by definition of G{z), the 
expansion of G{z) is given by 



G{Z) = J^anZ"", 9n = {%\{CgCpY\%). (3.36) 
n=0 

From (I3.32P and (13.360 . one immediately finds 

0+(O)=^?fc. (3.37) 
Thus, with fl3:35|) and ([SJZD, fl2:23|) and imply novel relations between 0^(g) and 0^(g) 



2n 



^^^(^) = 4jKm E(-l)''^+(^)^"' (3-38) 



' q\ 2n-l 



4>n<l) = . . E (-1) Vt(?)^72n-i-.- (3.39) 

4smn [^) 

Recall that 0+(q') and 0+(q') are related through the Fourier transformation originally. 
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4 Exact results 



In this section, we give some exact results for the grand partition function at k = 1. After 
the study of the algebraic treatment in the previous section, our strategy to compute Tr 
explained in Section |2] is improved. There are three steps to know Tr p". We first fix by 
using the recurrence relations (12.251) . We then compute Tr p" separately for odd n and for even 
n. The computation of Tr p2"+i reduces to the problem to compute the expansion coefficients 
Qn of G{z), which are fixed by </'+(g'). For the computation of Tr p^", on the other hand, we use 
the formula (I2.24p . From this prescription, we have determined the exact partition functions 
of the ABJM theory at A; = 1 up to iV = 9. 



4.1 Exact partition functions 

Let us start by determining 4>^{q). Since the integral equations (I2.25P have the convolution 
integral, the Fourier transformation is useful to solve (I2.25p . Let us consider 4>^{q). The initial 
condition is (f)^{q) = 1. One can trivially find 4>\{q) = Cg/2. Using the Fourier transformation, 
the convolution between Cq and (p^iq) can be evaluated as 

^^,-.^U.')4/|e*«cJ^±,.„ (4.1) 

where we have used the formula 

Cp = J (4-2) 

Thus we obtain 4>\{q) = gs2g/(47r). Iterating this method, we can fix 0'^(g) recursively. We 
list the explicit forms of the first seven functions below 

11 1 

</>+(g) = 1, </>+(?) = -Cg, </)+(g) = -^(lS2g, (pliq) = Yo^CginCg - 2qs2g), 

= J7^{'^q'4, + -^'{cg - - 4/2)), 



1287r2 

^+(?) = ^^Cg{4qhlg + A7TCgiqs2g + l)- 11^34 + 4/2)), 



+ ^('?^.-lK4 + I^c^(l + 2c.). (4.3) 
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Similarly (f)^{q) is also fixed. The explicit forms for < A; < 6 are given by 

0° (g) = 1, <PH<i) = = - 

-(^) = 4-^' (4 - 2gs,(l + 2c,) - 7r(l - 2c,)), 



a3 



levr 

-^^^ = 3^^''^. + - c.) - 2g.2,(l + 2c,)) - Yi^c,c2/2(6c, + 5), 

1 ' 2„2 2 2 I o„„3„2 , 2 2 o„2\ 



+ 7^-^9/2 (12CgC^/2 + 8c,C^/2 + 5Cg/2 ~ Sc^/g)' 



256 

'^-^^^ = 384^^''^''? + 128^^~^'''''''''/2 + + 4gs,c^sJ/2 + 4gs,c3 

+ 156c^c;/2 + 26c;c2/2 - 24c3 - 4c; + 456sjc^ - 120c^s;/2) 
+ 1^(288^^2 - 724 + 52qslcl - Ugs^cl - ISgs^c, + IGOsJcJ - 36sjc, + 20c;) 

+ Y^(-32s;c;c;/2 + 224s;c; - 54c;s;/2 + 72c3c;/2 + 12c;c;/2 - Qc^ - cj). (4.4) 

Now we proceed to computing Trp^""*"^. For this purpose, let us start by computing 
gn = (Og|(c,Cp)'^|Op). It is easy to see that gn has the following integral form 

9n = l ^(0,|(c,c,)-|g)(g|(c,c,)— 10,) = J ^~^'^{q)<Pr-{q). (4.5) 

Note that Qn obviously does not depend on the insertion position m (m = 0, ■ ■ ■ , n), thus we 
can choose it handily. If m = is chosen, we obtain f l3.37p because of (t)%{q) = 27r6{q). Using 
the results (14.31) . we find the first fourteen values of (?„ 

1 1 71-2 1 -7r2 + 37r + l 



^^-4' ^^"Stt' QAtt' ''~2^7r^' ~ 2%^ 

97r3 - 287r2 + 6 _ -4571^ + 927?^ + 1627r - 12 

211 ■ 3%3 ' ^7 - 214 . 32^3 ' 

IStt^ - 567r2 + 3 _ 727r4 - 2197r3 - 747r2 + 1627r + 3 

215 . 3%4 ' 39 - 217 . 3%4 ' 

-24757r5 + 54847?^ + 76507r3 - MOOtt^ + 30 



^11 

912 



219 ■ 32 ■ 5%^ 

49507r5 - 86847r4 - 227257r3 + 17007r2 + 60757r - 30 
222 . 32 . 52,^5 ' 
20257r6 - 571507r4 + 11831271^ + 1323007r3 - 840071^ + 90 



223 . 34 . 52^6 

_ -16807571^ + 4936687r5 + 2731127r4 - 52245071^ - 97507r2 + 656107r + 90 

- 225 . 34 . 5%6 • ^4-6) 
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These lead to the exact value of Tr p^n+i from f l3.3ip 

1 rr.^ 7r-3 10 -TT^ 

Trp = -, Trp = , Trp — 



4' IGtt ' 2567r2 

7 _ -277r3 + 497r2 + 1267r - 42 ^ 9 _ Svr^ - 207r2 - 9Qtt + 12 

47257r5 - 63037r^ - 253007r3 - 1210071^ + 231007r - 660 



Tr p 
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217 . 32 . 527^5 



13 _ -303757r6 + IOIMtt^ + 6552007r3 + 9789007r2 - 5616007r + 4680 
""P - 221 . 34 ■ 5%6 • ^ 

Finally let us see the computation of Trp^". We have already determined (f)±{q)- Using 

these results, we can know the matrix elements p±{qi, q2)- The explicit forms of p±{qi, q'2) are 

given in Appendix O These matrix elements enable us to compute TrpJ. Here we give the 

exact values of Trp^" (1 < n < 4) 

2 1 4 7^2-8 6 97r-28 

^^^+=16^' ^^^+=512^' ^^^+=2^^' 

Trp-^^I^^±1?^, (4.8) 



and 



TT 



217 . 37r2 



ll7r2 - 327r-8 _ « -Il7r2 + 47r + 96 



Tr p_ = — — — , Trp_ = , Trp' 



TrP- = ^^W-^s • (4-9) 



167r ' 5127r2 ' 212712 

-14677r3 + 35207r2 + 53767r - 6144 
217 . 3271 

Note that Tr p" is computed by using fl2.12p . 

Collecting the above results, the grand partition function is expanded as 

_ z ^2 TT - 3 3 10 - 7r2 4 -97r2 + 207r + 26 5 3671^ - 12l7r2 + 78 g 

+ ^ + ^ + +^T^^^ + ^2 ^ + 214 . 32^3 ^ 

-757r3 + 1937r2 + 1747r - 126 7 105371^ - 201671^ - 41487r2 + 876 § 

^ 216 . 37r3 ^ ^ 221 . 327r4 ^ 

55 1 77r4- 134807r3- 153487r2 + 88807r + 4140 „ .^n 

+ 223 ■ 3%4 ' + )• (4-10) 

Let us recall that the coefficient of z'^ in is just the partition function Z{N), thus we 
obtain the exact ABJM partition function Z{N) at /c = 1 up to TV = 9. The partition functions 
for higher TV can also be computed in this way. 



4.2 Novel relation between G{z) and H±(z) 

From the results in the previous subsection, we find a novel relation 

00 2n— 1 °° 2n 

logG(z) = 1^ Tr p2"-i + |;7(Trp+ " Trp"") + 0(^^°), (4.11) 

n=l n=l 
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in the expansion around z = 0. As we will see below, it is shown that the above relation holds 
exactly for all z, that is, 

G{z) = ^!Mi±^ = (4.12) 

Before proceeding to the proof, let us comment on a consequence led by the above equality. 
Using f l4.12p . the grand partition function is written as 

E(z) = det(l - z^pI)G{z). (4.13) 

Since the expansion coefficients of G{z) are just given by 0+(O) (see ( I3.37P ). the function 
G{z) is fixed by 0+(O). On the other hand, the matrix elements p^{qi, g'2) are fixed by 4>^{q) 
through f l2.24p . Therefore det(l — -z^p+) is also determined by (t>^{q). We conclude that the 
total grand partition function is completely fixed by (f>^{q) (or p+) only! We do not need 
any information about (t>'l{q) (or p_) to know ^(2;). The novel relation we have found in 
this section suggests that we have an alternative strategy for computing the grand partition 
function. Previously we need to compute Trp". In particular, to compute Trp^", we need 
both of Trpj." and sum up them: Tr p^" = Tr p^" + Trp^". After knowing the relation, the 
ingredients we need are the coefficients Qn of G{z) and Tr p^". Tr p^" is automatically fixed by 
the above relation. 

The rest of this subsection is devoted to the proof of f l4.12p by using the properties of the 
parity-preserving Hankel matrix M. Using the relation fl3.3ip . one can check that the odd 
power part of fl4.12p holds for all-order in z. Thus we concentrate ourselves to the even power 
part. The statement to be shown is now rephrased as 



det(l - z^pl) 
det(l - z^pD' 



G{z)G{-z) = (4.14) 



Let us recall that M decomposes into two Hankel matrices M±. From fl2.38p . we find 

Ml = M+T_T+M+ = M+PM+, (4.15) 

where P = 1 — |eo)(eo| and |eo)(eo| is the projector to the zero-th component. This relation 
implies 

TrM!" = Tr(A/2p)'^ (4.16) 

The right hand side looks similar to fl3.18p and therefore suggests that a similar manipulation 
is possible. All we have to do is to rewrite the right hand side by expanding P = 1 — |eo) (eo|. 
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Again, this is most easily done by studying the determinants. Let us rewrite the right 
hand side of (HTTil) . We use the fact ^±(2) = det(l + zM±), then 

^ ^ ' det 1 + -±^|eo)(eo| . (4.17) 



det(l-z2M2) det(l-^2M2) ^ l-z'^Ml 

From the similar argument to derive (13.231) . this determinant becomes simple 



Expanding this in z, we obtain 

det(l - z'^Ml\ 



det(l - z'^Ml 



X]^2"(M^«)oo. (4.19) 



+ ' n=0 

Now we want to write the right hand side of (14.191) in terms of the coefficients Qn of G{z). To 
do this, we use (12.471) and (I3.38p . Setting g = in the second equation of (I2.47p . we obtain 

40^(0) = (M^-i)oo. (4.20) 

On the other hand, if setting g = in (I3.38p . we also obtain 

2n 

40^-+i(O) = Y.{-irgr,,g2n-^. (4.21) 

m=0 

where we have used (13.371) . Comparing these relations, we find 

2n 

(Mf )oo = Y,{-irgmg2n-m. (4.22) 

m=0 

Substituting (I4.22p into (I4.19p . one can confirm that (I4.14p really holds. 



5 Numerical results 

In the previous section, we compute the (grand) partition function analytically. Here we 
discuss another approach to analyze the partition function of the ABJM theory. The key idea 
is to use a link between the Fredholm theory for the kernel like (12. 4p and TEA system. This 
novel link was first considered by Zamolodchikov to study two-dimensional polymers in [25] . 
The kernel appearing there is similar to (12.40 . Thus we can use some of the results in [25] . 
Particularly, we can write down the TBA-type integral equations that give Tr p". Though 
the TEA equations can not be solved analytically in general, it is easy to solve them at least 
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numerically. In this section, we numerically solve the TBA-like equations and compute Trp". 
We compare them with the exact results in the previous section. Our numerical results show 
the very good agreement with the exact ones, and further give a prediction of Tr (and also 
the partition function), of which we have not found the exact values, with a very high degree 
of accuracy. We also study numerically the eigenvalue spectrum of the density matrix. 

5.1 Partition function from TBA-like equations 

In order to see the relation to the TBA-like equations, it is convenient to rewrite the kernel 
( 12. 4p as the form 

g-|C/(9i)-|C/(g2) 

2cosh(^) ' ^'-'^ 

where 

U{q) = log 2 cosh I . (5.2) 
In [23], Zamolodchikov considered the following integral 

/oo roo „-2u{xi) 2u{xn) 

.oo J^oo cosh(^) cosh(2i^) . . . cosh(2ys) 

One can easily check that if u{x) = U{kx)/2 is chosen, the integral of Rn{x) is related to 
Tr p""*"^ thorough 

/oo 
rfxi?„(x) = (47r)"+iTrp"+^ (5.4) 
-oo 

Let us introduce the generating functional of Rn{x), 

/ Z \ 2m / ^ \ 2m+l 

Re{x\z) = ^2m(a;), Ro{x\z) = ^ ^2m+i(a;), (5.5) 

m=0 m=0 

R{x\z) = Re{x\z) \ Ro{x\z) . (5.6) 
From (12. 6 p and (15. 4p . we obtain {^z = e^) 

— / dxR{x\z) = -logS(^) = e"'^J'(p), (5.7) 

where J(p) = logS(2;) is the grand potential. Thus if we know the function R[x\z), we can 
obtain the grand partition function. 
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One of the important results in [25] is that the function Re{x\z) is conjectured to be the 
solution of the following TBA-type integral equations 

/oo 
dx' Ki{x — x') log(l + ri'^{x'\z)), (5.8) 
-oo 

POO 

r]{x\z) = —z / dx' Ki{x — x')Re{x'\z), (5.9) 



oo 



where 7]{x\z) is an auxiliary function, and we have introduced the following functions 

2 cosh I 27rcoshx 

These two equations determine the functions Re{x\z) and ri{x\z) uniquely for given z. Once 
Re{x\z) and ri{x\z) are determined, the function Ro{x\z) is given by 

/oo 
dx' K2{x — x') aj:ctanrj{x'\z), (5-11) 
-oo 



Mx) = — -7T-. (5.12) 



where 

1 

vr cosh^ X 

Therefore we can obtain the whole function R{x\z) from these integral equations. In our case, 
the driving term in (15.81) takes the form — log(2 cosh(/cx/2)), which is not a standard form 
— m cosh X in ordinary TBA systems for 2d relativistic quantum field theories. Note that this 
conjecture has been proved in 



5.2 Numerical results from TBA 

Now we want to solve the TBA-like equations (15.81) . (15.91) and (15.111) and to compute Trp". 
For this purpose, we expand ri{x\z) as 

r]{x\z) = [-^) V2m+i{x). (5.13) 

m=0 

Substituting (15.51) and (15.131) into (15. 8p . we obtain the following first few equations, 

logi?o(a;) = log s{kx), = Ki*r]f, 



Rn x) 2 \Rn x) / - I 2 
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where the convolution is defined by 



oo 



f *9 



dx' f{x — x')g{x'). 



(5.15) 



— oo 



Similarly, from (15. 9p . we obtain 



(5.16) 



for m = 0, 1, ■ ■ ■ . Starting with Ro{x) = s{kx), one can determine ?7i(x) from (I5.16p . then 
R2{x) is determined by the second equation in (I5.14p . The fiow of the determined functions 
is as follows, 



The integral equations (15.141) . (I5.16p . (I5.18P etc. are solved recursively. In principle, we 
can solve them analytically by evaluating the convolutions of the known functions at each 
step. Technically, however, it is very messy to do that. Alternatively, here we solve these 
sequential integral equations numerically at A; = ij^l If Rn{x) is determined, we can compute 
Trp"+^ from 

We have indeed computed the numerical values of Tr p" from such integral equations. In 
Tabled], we summarize the numerical values for 2 < n < 16. The numerical results show the 
very good agreement with the exact values obtained in the previous section. We would like to 
emphasize that the TBA-like equations are valid for general k (even for non-integer values). 
The only difference is the initial condition Ro{x) = s{kx). We can predict the high-accurate 
values of the partition functions for various and k in the same manner. 

In this subsection, we expand all the unknown functions Re{x\z), Ro{x\z) and ri{x\z) in z 
to compute Trp". We note that the TBA-like equations (15. Sp . (15. 9 p and (15. lip can be also 
solved for given z by iteration as we will see in the next section. 

5.3 Eigenvalue distribution 

Here we study the eigenvalue distribution of the density matrix numerically. For this purpose, 
we use the Hankel matrix representation discussed in Section [2^21 Since this matrix is discrete, 

■^We have also analytically checked that the TBA-like equations indeed reproduce the correct Tr p" up to 
n = 5. 



Ro{x) r]i{x) R2{x) 773(0;) Ri{x) • • • . 



(5.17) 



In the similar way, Rn{x) for odd n is fixed from (15. lip such as 




(5.18) 



25 



Table 1: Numerical values of (27r)"Trp" at = 1. These numerical values are computed by 
using the TBA-like equations (15. 8p . (15. 9 p and (15. lip . Our numerical results nicely agree with 
the exact values in the previous section. 



n (27r)"Trp" Errors 



2 





896604773477443035477301058329286341729730 


1.1 


X 


10" 


-40 


3 





698731738515872159486421033736470898141922 


1.7 


X 


10' 


-40 


4 





587330256906414478771564191224945948765499 


2.1 


X 


10' 


-40 


5 





505385252214593561502730078447293192743804 


2.2 
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10' 


-40 


6 





438346955317146787321292739552212309497233 


2.3 


X 


10" 


-40 


7 





381274926859615816279007202371073982670788 
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10" 


-40 


8 





331972050713896807160802466534622236608778 


2.2 


X 


10" 


-40 


9 





289152279521382289292363666349330742224258 
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X 


10" 


-40 


10 





251890109244736461833445267225791497711006 










11 





219440841867035365985490063842951367327288 


1.7 


X 


10' 


-40 


12 





191175332952772245351675601074802777111664 










13 





166551755978085254984697038651649345830601 


1.4 


X 


10' 


-40 


14 





145100086034249146678243847815416397520091 










15 





126411483244355445082928632866179620259652 










16 





110129973053738524406779285719767961758381 











we can evaluate its eigenvalues by the level truncation. Furthermore the Hankel matrix ele- 
ments Mij only depend on i + j. This means that we can reduce the 0(L^)-computation into 
(!?(L)-computation where L is the size of the matrix. This is the reason why we use the Hankel 
matrix M rather than the Chebyshev representation M'^ of the density matrix. Physically, 
we are interested in the energy spectrum of the Fermi gas system. Since the Hamiltonian is 
related to the density matrix as p = e~^, the n-th energy eigenvalue is given by 

En = -\0g\n, (5.19) 

where A„ is the eigenvalue of the density matrix. 

In Table |2l we summarize the numerical result of the energy eigenvalues of this quantum 
mechanical system, computed from the Hankel matrix M using the level truncation. In Figure 
[H we show the distribution of these energy eigenvalues. We have used the extrapolated values 
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Table 2: Numerical energy eigenvalues of the Fermi gas system at k = 1. The values at oo are 
obtained by extrapolating the data with the functions e*-'^-' + e^^^/L + e^^-*/-^^ + e'-^V-^'^ where 
L is the size of the matrix. E±^n mean the n-th energy eigenvalues of the systems described 
by the density matrices p±. 



Size 


Eq (-E+,o) 


El {E^,o) 


E2 


Es {E^,i) 


Ei {E+^2) 


^5 (^-,2) 


Eg (-E+,3) 


100 


1.975939 


2.977194 


3.739240 


4.432293 


5.142744 


5.895999 


6.709097 


200 


1.975808 


2.974644 


3.721872 


4.372444 


5.003785 


5.656839 


6.350546 


400 


1.975772 


2.973812 


3.715146 


4.344395 


4.928247 


5.510457 


6.116287 


1000 


1.975760 


2.973523 


3.712410 


4.330522 


4.883362 


5.408403 


5.933915 


2000 


1.975759 


2.973474 


3.711873 


4.327247 


4.870444 


5.372886 


5.860046 


4000 


1.975758 


2.973460 


3.711711 


4.326134 


4.865368 


5.356542 


5.820704 


10000 


1.975758 


2.973456 


3.711658 


4.325733 


4.863290 


5.348715 


5.798487 


oo 


1.975758 


2.973449 


3.711541 


4.324808 


4.859445 


5.339096 


5.782347 



as in Table [2]to plot this graph. We observe that the square of the energy eigenvalue £"„ shows 
the good linear behavior on the excitation level n for small n: 

El - El ^ yn. (5.20) 

Surprisingly, the coefficient of n in this equation agrees with the one appearing in the relation 
between the number of states and the energy in the thermodynamic limit [17]: n{E) = + 
no. 

Let us compare these results with the estimations from the exact results in the previous 
section. We can estimate the ground state and the first excited state energies from Tr 
because 

00 

Trp^ = $^A^^,,^A^,o, (5.21) 
for large n. Using (14.81) and (lC.7p . we find 

Eo = ^+,0 ~ - log = 1.9757623 ■ ■ ■ , (5.22) 

which matches the estimation from the Hankel matrix in Table [2l The order of the subleading 
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n 

2 3 4 5 6 



Figure 1: The (square of) energy spectrum of the Fermi gas system against the excitation 
level n. The dots correspond to the extrapolated numerical eigenvalues in Table [2] while the 
sohd line is the plot of E'^ — Eq = 7r^n/2. 

corrections is 0{\l/\l) ~ (9(10~^). Similarly, 

E, = E_,o ~ - log [t^J = 2.9735127 ■ ■ ■ , (5.23) 



where we have used the exact results (14.91) and fIC.Sp . The result is again consistent with the 
value in Table |2j The subleading corrections come from 0{\l/\l) ~ C(10^^). 



6 Comparison with perturbative results 

In this section, we compare our results with the perturbative ones in |17] . 

We first consider the free energy of the ABJM theory. In [13 HZ], it was shown that the 
ABJM partition function is written in terms of the Airy function if one neglects the instanton 
corrections. Since we have obtained the exact partition function Z{N) up to = 9 in Section 
m we can compare them with the perturbative results. This comparison allows us to know the 
behavior of the non-perturbative effects. Our results suggest that the leading non-perturbative 
correction scales as e"^'^^^^. 

We next consider the expectation value of A^ (or the derivative of the grand potential). 
Using the TBA-like equations in the previous section, we can numerically compute (A^) as a 
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function of the chemical potential /i. Our numerical results show the good agreement with 
the perturbative solution in [T7] in the large fi limit. 



6.1 Free energy 

Let us consider the free energy of the ABJM theory. As was shown in [151 [13 , the partition 
function is written in terms of the Airy function if one neglects the instanton corrections, 

Z,iN) = zi''"^\N) + 4°P)(iV), (6.1) 



-1/3 



'^/^ / 1 k 
TT^k V 3k 24 



(6.2) 



where Z'j^'^\N) is the non-perturbative contribution. The function A{k) has the following 
small k expansion, 

In [31], it was conjectured that this function has the following integral expression, which should 
be valid for arbitrary k, 

A(k) = - log 2 - _ 1 log A; + - log f-) + 2C'(-1) 

3 Jo e''^ — 1 X xsinh'^x^ 



7 {-^ 7^ • (6-4) 



Figure [2] (a) shows the behavior of the free energy F{N) = — log Z{N) at A; = 1. The result 
(16. 2 p shows the very close behavior of the exact results even though the non-perturbative 
corrections are dropped. It is surprising that these corrections at /c = 1 is very small even for 
small A^. 

Let us consider the non-perturbative corrections to the free energy 

^{-P) (AT) = t) (AT) _ ir{A-y) ^ ^g_5) 

where fI^^"^\n) = - log zI^^^\n)^ We find ^^^"^^(A^) is negative for all A^ < 9. In [H [JT] 
it was noticed that there are two kinds of non-perturbative instanton effects contributing to 
the free energy. One of them is the worldsheet instanton exp(— 27rA/2A^/A;), while the other 
is the membrane instanton exp(— 7rv2A?iV) which comes from D2-branes wrapping on MP^ in 
CP . Since in the case k = 1 there is no RP^ to wrap, we expect that the membrane instanton 
effect vanishes and the non-perturbative corrections start from exp(— 27r-\/2iV). 



^In the definition here, J^i"''^(iV) 7^ - log 4"^' (TV) where = " zi^"^\N). 
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N 



(a) (b) 

Figure 2: The behavior of the free energy at k = 1. (a) The dots show our exact values while 
the solid line shows the free energy computed from (16.21) . (b) The (log of) non-perturbative 
correction Fi^'^\n). The dots are the exact data, and the solid line is the fitted function (16. 6p 

with f ieri) . 

In Figure [2] (b), we tried to fit our exact data of log \ Fi^^\n)\ by the function 

log iFi"""^ {N)\ a + bhgN - cy/N, (6.6) 

and found the best values of the fitting 

a = -3.17, 6 = 2.01, c = 9.00 . (6.7) 

Note especially that the value of c is very close to 2n\/2 = 8.88577. This result strongly 
suggests that the membrane instanton e"'^^^'^^ is absent when k = I, and the leading non- 
perturbative correction behaves as F^^^N) ~ ^-"^-^V^n large A^. 

Finally, let us briefly comment on the pre-factor of the instanton correction e"^'^^'^. From 
the general argument in [17], one expects that F^^^\n) scales as 

(AT) ^ (AN + bVn + C)e^27rV2]v ^q q^ 

with some coefficients A, B and C . Our analysis ( 16. 7p seems to indicate that the leading term 
of F^^^\n) behaves as jv^e"^'"^^ rather than Ne~'^'^^^. However, our result up to = 9 
is insufficient to tell. One possibility is that ^^^'^^^(A^) may actually scale as (16. 8p for much 
larger N than = 9. It would be interesting to study the behavior of the non-perturbative 
correction Fi^'^\n) in more detail. 
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6.2 Expectation value of 

In the previous section, we solved the TBA-hke equations by expanding all the unknown 
functions in z in order to compute Tr p". It is possible to solve them for given z by iteration. 
Let us compare the results from such integral equations with the perturbative ones as functions 
of 2; = e'^. For this purpose, it is convenient to focus on the derivative of the grand potential 
J(yu) because this quantity naturally appears in the TBA approach as in (15. 7p . From the 
standard argument of statistical mechanics, the derivative J'(/i) corresponds to the expectation 
value of in the grand canonical ensemble 

{N) = J'ifi) = — / dxRixle"), (6.9) 

where (15. 7p has been used in the last equality. Using (12. 4p . the grand potential is also expressed 
in terms of Tr p". The expectation value (A^) is thus written as the expansion around p = — 00 

00 

(AT) = -^(-e^)"Trp". (6.10) 

n=l 

Let us consider the perturbative contribution of J(p). As shown in [17], the grand potential 
of the ABJM theory is given by 

if non-perturbative corrections are neglected. Thus the derivative of J{fi) is given by 

Note that, in (I6.12p . the contribution from A[k) is dropped. 

We numerically solve the TBA-like equations (15. Sp . (15. 9 p and (15. lip at k = 1 for various 
values of 2; = e^, and compute the expectation value of A^ by using (16. 9p as a function of p. 
The results are shown in Figure O In the figure, the dots represent our numerical results from 
the TBA-like equations. The red solid line is the plot of the truncated expansion (I6.10p up 
to n = 9 where we have used the exact values of Tr p". The green dashed line represents the 
perturbative result (I6.12p . We find that our numerical results are close to the perturbative 
result (I6.12P in the large p limit while in the limit p — — 00 the expansion (I6.10p agrees with 
the numerical results as expected. The TBA interpolates both results smoothly. This is not 
surprising because the TBA results contain all the non-perturbative corrections, which are 
neglected in (I6.12p . 
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Figure 3: The behavior of the expectation value of as a function of /i. The dots show the 
results obtained from TBA-like equations. The red solid line represents the finite-truncated 
expansion f lG.lOp up to = 9 with the exact Tr p" while the green dashed line represents the 
perturbative result f l6.12p . 

Finally, let us give a comment on an underlying technical issue in this TBA system. As 
was mentioned in [25], the iterative solution of the TBA equations describing the Fredholm 
determinant converges only when z is in the regime \z\ < I/Aq where Aq is the largest eigenvalue 
of the kernel p. li \z\ crosses over this value, the iteration method does not work any more. In 
our interested case k = 1, the bound of the convergence exists at \z\ = I/Aq ~ 7.21 (/x ~ 1.97). 

7 Discussion 

In this paper, we have studied the Fermi gas quantum mechanics associated to the ABJM 
matrix model. We have computed Tr p", which is the ingredient of the grand partition function, 
and obtained the partition function Z{N) up to = 9 as a result. In reducing the number 
of integrations by commuting coordinates and momenta, we have found an exact relation 
concerning the grand partition function, which is interesting by itself and very helpful for 
determining the partition function. We have also performed various numerical analysis. All 
of these results are consistent with our exact results and also with the perturbative analysis 
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in [I?]. 

Interestingly, all the partition functions we obtained are written in terms of polynomials of 
7r~^ with rational coefficients. Similarly, the functions (f)^{q) which lead to the determination 
of Tr p" are also expressed in terms of polynomials of tt"^, q, Sg, Cg, Sq/2 and Cg/2 with rational 
coefficients. Since we have computed the partition function exactly and integer numbers (or 
rational numbers) usually counts some information of the moduli space, we would like to know 
what kind of the information we are counting right now. 

We would like to understand better the origin of the anomaly in commuting the coordinate 
q and the momentum p. From the analogy of string field theory, typically these kinds of 
anomaly appear because of the singularity of the moduli space we are considering and give 
some information to the physical problem. For example, the associativity anomaly [32] in 
cubic string field theory [33] comes from the singularity at the boundary of the moduli space 
and relates to the translational invariance in string field theory. Also, the twist anomaly [33] 
in vacuum string field theory [3S| originates from the divergence at the midpoint interaction 
point and contains the information of the tachyon mass and the brane tension. 

We have shown that finding the energy spectrum of our Fermi gas system boils down to the 
diagonalization of the Hankel matrix M defined by fl2.39l) . As in the case of the diagonalization 
of Neumann matrices in open string field theory [3B] , the first step would be to find a (hidden) 
symmetry K commuting with M: [M, K] = 0. It would be very interesting to find the exact 
spectrum of M and write down the exact grand partition function. 

We showed in (12.131) that the grand partition function of the ABJM theory is factorized 
into two parts. A similar structure has appeared in the study of holographic dual of BPS black 
holes [37], and it was argued that going to the fixed charge sector leads to the entanglement 
of two sectors. Such entangled state has an interpretation of branched tree of baby universes 
and the Catalan number appears as the number of branchings [37]. It would be interesting to 
see if the factorized grand partition function of the ABJM theory has some interpretation on 
the bulk gravity side. 

We also would like to remark on an analogy with the so-called ODE/IM correspondence 
[551 In [301 1311 132], it was found that spectral determinants for certain Schrodinger 

equations satisfy functional relations. In [3H1EH], these spectral determinants were identified 
with Baxter's Q-functions of the corresponding integrable system, and the functional relations 
of the spectral determinants just correspond to the quantum Wronskian relations [33]. Our 
analysis implies that there is a non-trivial relation between the two determinants 'E± [z) because 
the total determinant S(z) is fixed by the information of either of the two sectors p±. It 
would be interesting to explore whether '^±{z) satisfy a functional relation like the quantum 



33 



Wronskian relation or not. 

Recently, the expectation values of some BPS Wilson loops of the ABJM theory were 
analyzed using the Fermi gas formalism [12] . It was found that the semiclassical expansion of 
such Wilson loop expectation values is divergent for k = 1,2. We hope that our formalism 
will give some improvements for the computation of Wilson loops. 
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A Recurrence relations for 

Let us derive the recurrence relations f l2.25p . Using f l2.22p . we find 

E±[q) J 2iT 2-K 
1 fdq' 



E±{q) J 271 

Using the identity 

1 



l^p^{q,q')E^{q')<Pi-\q'). (A.l) 



cosh q + cosh q' 2 cosh | cosh |- 
we find that the kernel is written as 

EM^'p^{q,q')E^{q') = '''^^ = c /^'^' + ^^'^^ 

2 (cosh q + cosh q') 2 

,22: 



riCg.g' + Cq + g>), (A.2) 



E-{qrP-{q,q')E-{q') = , , , ""^ \ = c,^^^^^^^^. (A.3) 



2 cosh ^ (cosh q + cosh q' 



Since 0j.(g)'s are even functions of q, we get 



<P'-{q) = I ^c,^i^^l±^tJ0^^(g') = c, I ^Cg^,tl4t-\q% (A.4) 
which are the recurrence relations ( 12.25^ . 
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B Hankel representations of (p^^q) and 

Let us derive fl2.47p . We expand and in tg = tanh(g/2) 



oo 

Cn 



= ^clY.tfal,'\ 4>i{q) = f (B.l) 

n=0 n=0 

We would like to show 

oo 
m=0 

Since 0° (g) = 1 and (t>\{q) = Cq/2, one immediately finds 

af = l, /3« = 5.o. (B.3) 
Let us substitute (IB.ip into the recurrence relation fl3.33p . Using the expansion 

oo 

c,_,, = 2c,c,, J]t^t^„ (B.4) 

n=0 

we obtain 

0t (?) = / 1^ f 2^.^.' E ■ 44 E tl^af^ (B.5) 
\ n=0 / e=o 

The integral over q' causes the Hankel matrix because of fl2.29p . Thus we find 

00 

</>+(?) = 4c; EC(^+)"^"?"'^- (B-6) 

e,n=o 

Comparing this equation with (IB.ip . we obtain 

00 

«W = ^(M+)„,af-^). (B.7) 
e=o 

From this relation and the initial condition (IB.3p . one can show that an^ is given by ( IB.2p . 
The derivation of (3n^ is straightforward. 
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C Explicit forms of ^±(^1,^2) 

Let us summarize the explicit forms of p±{qi, ^2) up to n = 4: 

2 1 1 1 

P+(gi, 92) = ^cosh(f ) cosl(fy cosh(f ) + cosh(f ) ' ^^-^^ 

^3. 1 V/C0sh(f)c0sh(f) f 02 A .(.21 



P+(9l,92) 



IGtt cosh gi + cosh g2 \sinhgi sinhg2 2 cosh(Y) cosh(Y' 



1287r ^cosh(f )cosh(f) cosh(f) - cosh(f ) Vsinhgi sinhg2 



2 
1 



256cosh3/2(^)cosh3/2(f)^ 



and 



, 1 sinh(f ) sinh(f' 



167r ^cosh(^) cosh(^)(cosh qi + coshg2) 

qi . q2 TT 



sinhgi sinhg2 2 cosh(Y) cosh(Y' 
tanh(f )tanh(f ) 



pi (91,92) 



64cosh3/'(f)cosh3/2(f; 

1 tanh(f )tanh(f) / 1 tt 



2567r ^cosh(f ) cosh(f ) I sinh'(f ) sinh'(f ) cosh(f ) cosh(f 



1 tanh(f )tanh(f ; 



5127r ^cosh(f )cosh(f)(cosh(f ) - cosh(f )) 
gicoth(f) g2Coth(f) 



(C.3) 



2. 1 tanh(^)tanh(f) 

^"^^''^'^ 16 v/cosh(f ) cosh(f )(cosh(f ) + cosh(f )) ' ^ ' ' 



(C.5) 



cosh(f )sinh2(f ) cosh(f )sinh^(f ^ ' ■ ^^'^^ 

These matrix elements are directly derived from f l4.3p . f l4.4p by using fl2.23p and fl2.24p PI Using 
these matrix elements, we can also compute Tr Pi""*"^ up to n = 3 in addition to Tr p^" as given 
in (14. 8p and f l4.9p . The results are given by 

^ 3 _ 3-2^2 ^ 5 _ (5-8v^)7r + 20 



Trp+ = ^, Trp+ = , Trp+ 



64 ' ^+ 210 



vr 



(16v^- 21)7r2 - 147r + 28 



2i%2 



^ Similarly we also obtain p±{qi,q2) for n = 5,6,7. Since it can be obtained without difficulty, we shall 
omit their expressions here due to the complexity. 
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and 

2-^2 ^ 3 (l + 2v^)7r-12 ^ 5 (8^2 - 9)7r2 - 207r + 40 

^^^- = ^' ^'P- = eTn ' ^^^- = ' 

^ 7 -9(16^2 + 3)7r3 + 5187r2 + 7567r- 336 ,^ „, 

P- = ^^^T^, • (C.8) 

Interestingly, in the above results the terms with \/2 appear. Such terms are canceled out in 
Tr p^'"+^, and Tr p^n+i ^^ot contain them. 
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